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In this paper, we study modulus of continuity and rate of convergence of series of conditionally 
sub-Gaussian random fields. This framework includes both classical series representations of 
Gaussian fields and LePage series representations of stable fields. We enlighten their anisotropic 
properties by using an adapted quasi-metric instead of the classical Euclidean norm. We specify 
our assumptions in the case of shot noise series where arrival times of a Poisson process are 
involved. This allows us to state unified results for harmonizable (multi)operator scaling sta¬ 
ble random fields through their LePage series representation, as well as to study sample path 
properties of their multistable analogous. 
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1. Introduction 

In recent years, lots of new random fields have been defined to propose new models for 
rough real data. To cite a few of them, let us mention the (multi) fractional Brownian 
fields (see, e.g., [6]), the linear and harmonizable (multi)fractional stable processes [12, 36] 
and some anisotropic fields such as the (multi)fractional Brownian and stable sheets [3, 4] 
and the (multi)operator scaling Gaussian and stable fields [8, 9]. In the Gaussian setting, 
sample path regularity relies on mean square regularity. To study finer properties such as 
modulus of continuity, a powerful technique consists in considering a representation of the 
field as a series of random fields, using for instance Karhunen Loeve decomposition (see 
[1], Ghapter 3), Fourier or wavelet series (as in [5, 17]). This also allows generalizations to 
non-Gaussian framework using for instance LePage series [24, 25] for stable distributions 
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(see, e.g., [20]). Actually, following previous works of LePage [24] and Marcus and Pisier 
[29], Kono and Maejima proved in [21] that, for a G (0,2), an isotropic complex-valued 
Qf-stable random variable may be represented as a convergent shot noise series of the 
form 

+ 00 

n—1 

with {Tn)n>i the sequence of arrival times of a Poisson process of intensity 1, and (A„)„>i 
a sequence of independent identically distributed (i.i.d.) isotropic complex-valued random 
variables, which is assumed to be independent of (T„)„>i and such that E(jAij“) < 
-boo. When A„ = Vngn with (g„)n>i a sequence of i.i.d. Gaussian random variables 
independent of (14,, 7n)n>i, the series may be considered as a conditional Gaussian series. 
This is one of the main argument used in [7, 8, 12, 20] to study the sample path regularity 
of some stable random fields. Another classical representation consists in choosing Xn = 
Vn£n with (en)n>i & Sequence of i.i.d. Rademacher random variables that is, such that 
P(e:„ = 1) = P(en = —1) = 1/2. Both and e„ are sub-Gaussian random variables. Sub- 
Gaussian random variables have first been introduced in [17] for the study of random 
Fourier series. Their main property is that their tail distributions behave like the Gaussian 
ones and then sample path properties of sub-Gaussian fields may be set as for Gaussian 
ones (see Theorem 12.16 of [23], e.g.). In particular, they also rely on their mean square 
regularity. 

In this paper, we study the sample path regularity of the complex-valued series of 
conditionally sub-Gaussian fields defined as 

+ CO 

S{x) = '^Wn{x)gn, for x G Aid C (2) 

n—1 

with {gn)n>i 3- sequence of independent symmetric sub-Gaussian complex random vari¬ 
ables, which is assumed independent of (114)„>i- In this setting, we give sufficient as¬ 
sumptions on the sequence (II4)n>i to get an upper bound of the modulus of continuity 
of S as well as a uniform rate of convergence. Then, we focus on shot noises series 

+ 00 

S{a,u) = '^T~^^°‘Vn{a,u)gn, x={a,u) G A4+i C (0,2) x 

n—1 

with (Tn)^^i the sequence of arrival times of a Poisson process. Assuming the indepen¬ 
dence of iTn)n>ij (14)„>i and (ffra)„>i, we state some more convenient conditions based 
on moments of 14 to ensure that the main assumptions of this paper are fulfilled. In 
particular when Vn(a,u) := A„ is a symmetric random variable, one of our main result 
gives a uniform rate of convergence of the shot noise series (1) in a on any compact 
Ki = [a,b] C (0,2), which improves the results obtained in [II] on the convergence of 
such series. In the framework of LePage random series, which are particular examples of 
shot noise series, we also establish that to improve the upper bound of the modulus of 
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continuity of S, one has the opportunity to use an other series representation of S. On 
the one hand, our framework allows to include in a general setting some sample path reg¬ 
ularity results already obtained in [7, 8] for harmonizable (multi)operator scaling stable 
random fields. On the other hand, considering a as a function of it G we also inves¬ 
tigate sample path properties of multistable random fields that have been introduced in 
[13]. To illustrate our results, we focus on harmonizable random fields. 

The paper falls into the following parts. In Section 2, we recall definition and properties 
of sub-Gaussian random variables and state our first assumption needed to ensure that 
the random field S is well-defined by (2). We also introduce a notion of anisotropic 
local regularity, which is obtained by replacing the isotropic Euclidean norm of by 
a quasi-metric that can reveal the anisotropy of the random fields. Section 3 is devoted 
to our main results concerning both local modulus of continuity of the random held 
S dehned by the series (2) and rate of convergence of this series. Section 4 deals with 
the particular setting of shot noise series, the case of LePage series being treated in 
Section 4.3. Then Section 5 is devoted to the study of the sample path regularity of 
stable or even multistable random helds. Technical proofs are postponed to Appendix 
for reader convenience. 


2. Preliminaries 

2.1. Sub-Gaussian random variables 

Real-valued sub-Gaussian random variables have been dehned by [17]. The structure of 
the class of these random variables and some conditions for continuity of real-valued 
sub-Gaussian random helds have been studied in [10]. In this paper, we focus on con¬ 
ditionally complex-valued sub-Gaussian random helds, where a complex sub-Gaussian 
random variable is dehned as follows. 

Definition 2.1. A complex-valued random variable Z is sub-Gaussian if there exists 
s G [0, -boo) such that 

VzGC, (3) 

Remark 2.1. This dehnition coincides also with complex sub-Gaussian random vari¬ 
ables as dehned in [15] in the more general setting of random variables with values in 
a Banach space. Moreover, for a real-valued random variable Z, it also coincides with 
the dehnition in [17]. Kahane [17] called the smallest s such that (3) holds the Gaussian 
shift of the sub-Gaussian variable Z. In this paper, if (3) is fulhlled, we say that Z is 
sub-Gaussian with parameter s. 

Remark 2.2. A complex-valued random variable Z is sub-Gaussian if and only if Ift(Z) 
and Q(Z) are real sub-Gaussian random variables. Note that if Z is sub-Gaussian with 
parameter s then E(5R(Z)) = E(7y(Z)) = 0 and E(5R(2')^) < as well as E(Sf(Z)^) < s^. 
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The main property of sub-Gaussian random variables is that their tail distributions 
decrease exponentially as the Gaussian ones (see Lemma A.l). Moreover, considering 
convergent series of independent symmetric sub-Gaussian random variables, a uniform 
rate of decrease is available and the limit remains a sub-Gaussian random variable. This 
result, stated below, is one of the main tool we use to study sample path properties of 
conditionally sub-Gaussian random fields. 

Proposition 2.1. Let (3n)„>i be a sequence of independent symmetric sub-Gaussian 
random variables with parameter s = 1. Let us consider a complex-valued sequence a = 
(an)n>i such that 

+ 00 

||a||^2 = ^ \anf < -1-00. 

n—1 

1. Then, for any t £ (0,-foo), P(supjvgN\{o} II]n=i “"SnI > t\\a\\p) < 8e“‘ /®. 

2. Moreover, the series converges almost surely, and its limit Jf,n=i^rign is a 

sub-Gaussian random variable with parameter ||a||^ 2 . 

Proof. See Appendix A. □ 

Remark 2.3. In the previous proposition, assuming that the parameter s = 1 is not 
restrictive since a„ can be replaced by a„s„ and pn by gnjsn when is sub-Gaussian 
with parameter s„ > 0. 

2.2. Conditionally sub-Gaussian series 

In the whole paper, for d > 1, it'd = 11^=1 C is a compact d-dimensional interval 

and for each integer N €N, we consider 

N 

SNix) = '^Wn{x)gn, X&Kd, (4) 

n—1 

where = 0 by convention and where the sequence (Wn, 5n)n>i satisfies the following 

assumption. 

Assumption 1. Let {gn)n>i xnd (W,i)„>i be independent sequences of random vari¬ 
ables. 

1. ign)n>i is a sequence of independent symmetric complex-valued sub-Gaussian ran¬ 
dom variables with parameter s = I. 

2. {Wn)n>i is ® sequence of complex-valued continuous random fields defined on Kd 
and such that 

- 1-00 

almost surely E \Wn{x)\'^ < -1-00. 

n—1 


'ix G Kd, 
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Under Assumption 1, conditionally on each Sat is a sub-Gaussian random 

field defined on Kd. Moreover, for each x, Proposition 2.1 and Fubini theorem lead to 
the almost sure convergence of Sn{x) duS N ^ +oo. The limit field S defined by 


+ 00 

S{x) = '^Wn{x)gn, xGA^dCR'^, (5) 

n—1 


is then a conditionally sub-Gaussian random field. In the sequel, we study almost sure 
uniform convergence and rate of uniform convergence of as well as the sample 

path properties of S. 

Assume first that each gn is a Gaussian random variable and that each Wn is a deter¬ 
ministic random field, which implies that 5" is a Gaussian centered random field. Then, 
it is well known that its sample path properties are given by the behavior of 


s{x,y) 


' +00 


1/2 




\n—l 


x,y£ Kd, 


( 6 ) 


since is proportional to the variogram {x,y) ^ v{x,y) := E[|5'(x) — S'(?/)|^]. In the 
following, we see that under Assumption 1, the behavior of S is still linked with the 
behavior of the parameter s. In this more general framework, a key tool is to remark 
that conditionally on S' is a sub-Gaussian random field and the random variable 

S(x) — S{y) is sub-Gaussian with parameter s{x,y). 

We are particularly interested in anisotropic random fields S (and then anisotropic 
parameters s). Therefore, next section deals with an anisotropic generalization of the 
classical Holder regularity, that is, with a notion of regularity which takes into account 
the anisotropy of the fields under study. 


2.3. Anisotropic local regularity 

Let us first recall the notion of quasi-metric (see, e.g., [32]), which is more adapted to 
our framework. 

Definition 2.2. A continuous function x —>■ [0,-|-oo) is called a quasi-metric 

on if 

1. p is faithful, that is, p(x, y) =0 iff x = y; 

2. p is symmetric, that is, p{x,y) = p{y,x); 

3. p satisfies a quasi-triangle inequality: there exists a constant k> 1 such that 

Vx, y,zG p{x, z) < k(p(x, y) -f p{y, z)). 

Observe that a continuous function p is a metric on if and only if p is a quasi-metric 
on which satisfies assertion 3 with k= 1. In particular, the Euclidean distance is an 
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isotropic quasi-metric and its following anisotropic generalization 


{x,y)^ p{x,y) 


1/p 






where p> 0 and ai,... ,ad > 0, 


is also a quasi-metric. Such quasi-metrics are particular cases of the following general 
example. 


Example 2.1. Let us consider E a real d x d matrix whose eigenvalues have positive 
real parts and te ^ a continuous even function such that 

(i) for all x^O, te{x) > 0; 

(ii) for all r > 0 and all a; G TE{r^x) = rxEix) with = exp ((Inr)i?). 

The classical example of such a function is the radial part of polar coordinates with 
respect to E introduced in Chapter 6 of [31]. Other examples have been given in [9]. 

Let us consider the continuous function pE, defined on x R'^ by 

PE{x,y) = TE{x-y). 

Then, by definition of te, Pe is faithful and symmetric. Moreover, by Lemma 2.2 of [9], 
Pe also satisfies a quasi-triangle inequality. Hence, pe is a quasi-metric on R*^ and it is 
adapted to study operator scaling random fields (see, [7, 9], e.g.). 

Let us remark that since p^ defines a quasi-metric for E/j5 whatever /3 > 0 is, we may 
restrict our study to matrix E whose eigenvalues have real parts greater than one. Then, 
by Proposition 3.5 of [8], there exist 0 < H_< H <1 and two constants C 2 ,i,C 2,2 G (0,oo) 
such that for all a;, y G R'^, 

C 2 ,i min(||a; - yjj^, jja; - yjj^) < pe{x, y) < C 2 , 2 max(||x - yjj^, ||x - yjj-^), 

where || • || is the Euclidean norm on R"^. In [7, 8], this comparison is one of the main 
tool in the study of the regularity of some stable anisotropic random fields. Therefore, 
throughout the paper, we consider a quasi-metric p such that there exist 0 < El < H < 1 
and two constants C 2 ,i,C 2,2 G (0, oo) such that for all x,y G R'^, with jjx — yjj < 1, 

C2,i||a: - y||^ < p(x,y) < C2,2||a; - y||-- (7) 

Before we introduce the anisotropic regularity used in the following, let us briefly com¬ 
ment this assumption. 


Remark 2.4- 

1. The upper bound is needed in the sequel to construct a particular 2“^ net for p, 
whose cardinality can be estimated using the lower bound. 
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2. Using the quasi-triangle inequality satisfied by p and its continuity, one deduces 

from (7) that for any non-empty compact set Kd = 0^=1 U there exist 

two finite positive constants C 2 ,i{Kd) and C 2 ^ 2 (,Kd) such that for all x,y G Kd, 

C 2 .i{Kd)\\x - y\\^ < p{x,y) < C 2 , 2 {Kd)\\x - y\\^. (8) 

3. It is not restrictive to assume that H <1 since for any c > 0, is also a quasi-metric. 

We will consider the following anisotropic local and uniform regularity property. 

Definition 2.3. Let P G (0,1] and p G M. Let Xg G Kd with Kd C R^^. A real-valued 
function f defined on Kd belongs to 77p,ic^(xo,/3,p) if there exist 7 G (0,1) and C G 
(0, -boo) such that 

\fix) - f{y)\< Cp{x,yf\log{p{x,y))\'’ 

for all x,y G B(xo,j) CiKd = {z G Kd; ||z —xojj < 7 }. Moreover f belongs to 'Hp{Kd,P,p) 
if there exists C G (0, -boo) such that 

yx,yGKd, \fix)-f{y)\ <Cp{x,y)^[\og{l-Gp{x,y)-^)]"’. 


Remark 2.5. 

1. If / € 'Hp^K.i{xg,P,p), then / is continuous at xg. Moreover, since h^[log(l -b 
h~^)Y^ ~/t-s-o+ ^^|log(li)r and since p satisfies equation (7), / G Rp^Kdi^o,P,v) if 
and only if for some 7 > 0, f G 'Hp{B{xg,j) fl Kd,P, p). 

2. If / G 'Hp{Kd, P,p), then / G Rp^Kdixg, P,p) for all xg G Kd- The converse is also 
true since Kd is a compact. This follows from the Lebesgue’s number lemma and the 
boundedness of the continuous function / on the compact set Kd (see Lemma B.2 
stated in the Appendix for an idea of the proof). 

3. A function in Rp^Kd, /?, 0) may be view as a Lipschitz function on an homogeneous 
space [28]. Note also that when p is the Euclidean distance, for any /3 < 1 and p <0, 
the set Rp{Kd,P, p) (resp., Rp,Kd {xg,P, p)) is included in the set of Holder functions 
of order /3 on Kd (resp., around xg). 

4. Assuming /? < 1 is not restrictive since, for any c > 0, is also a quasi-metric. 

The introduction of the logarithmic term appears naturally when considering Gaus¬ 
sian random fields. Actually, [ 6 ] proves that for all /3 G (0,1], a large class of elliptic 
Gaussian random fields Xg, including the famous fractional Brownian fields, belongs 
a.s. to Rp^Kd{xg,P,l/2) with p the Euclidean distance (see Theorem 1.3 in [ 6 ]). More¬ 
over, Xiao [38] also gives some anisotropic examples of Gaussian fields belonging a.s. to 
Rp^Kdi^o, 1/2) for some anisotropic quasi-distance p = pe associated with E a diagonal 
matrix (see Theorem 4.2 of [37]). Finally, in [ 8 ], we construct stable and Gaussian random 
fields belonging a.s. to Rp^ Kdi^oA ~£j0) for some convenient px^ (see Theorem 4.6 in 
[ 8 ]). 
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3. Main results on conditionally sub-Gaussian series 

3.1. Local modulus of continuity 

In this section, we first give an upper bound of the local modulus of continuity of S 
defined by (4) under the following local assumption on the conditional parameter ( 6 ). 

Assumption 2. Let Xq & Kd with Kd = consider p a quasi¬ 

metric on satisfying equation (7). Assume that there exist an almost sure event Ll' 
and some random variables 7 > 0, /3 G (0,1], ij gM. and C G (0, + 00 ) sueh that on Lt' 

Vx,y G B{xo,-f)nKd, s{x,y) < Cp{x,y)>^\\og{p{x,y))\'^, 

where we recall that the conditional parameter s is given by (6). 

Note that the event fl', the random variables 7 ,/?, 77 , C and the quasi-metric p may 
depend on xq. 

Let us now state the main result of this section on the modulus of continuity. The 
main difference with [7, 8 , 21] is that we do not only consider the limit random field S 
but obtain a uniform upper bound in N for the modulus of continuity of Sn- 

Theorem 3.1. Assume that Assumptions 1 and 2 are fulfilled. Then, almost surely, 
there exist 7 * G (0, 7 ) and C G (0, -foo) such that for all x,y G B{xq,^*) n Kd, 

sup |S'jv(a:) - SN{,y)\ < Cp{x,yY\\ogp{x,y)\'^"^^^'^. 

Nm 

Moreover, almost surely {SN)NeN converges uniformly on B(xo,j*)nKd to S and the 
limit S belongs to Tip^Kdi^cn A- 1/2). In particular, almost surely S is continuous at 
Xo- 

Proof. See Appendix B.l. □ 

Strengthening Assumption 2, the uniform convergence and the upper bound for the 
modulus of continuity are obtained on deterministic set. Next corollary is obtained using 
some covering argument. 

Corollary 3.2. Assume that Assumption 1 is fulfilled. 

1. Assume that Assumption 2 holds for any xq G Kd with the same almost sure event 
LI', the same random variables fi and 77 , and the same quasi-metric p. Then Theo¬ 
rem 3.1 holds replacing B{xo,"f*) r\Kd by all the set Kd and almost surely S belongs 
to np{Kd,fi,p + i/2). 

2. Assume now that Assumption 2 holds with a deterministic 7 . Then Theorem 3.1 
holds replacing B{xo,j*)nKd by B{xo,j)r]Kd and almost surely S belongs to 
TLp{B{xo,j)r\Kd,fi,ri +1/2). 
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Proof. See Appendix B.l. □ 

When considering S an operator scaling Gaussian random field, note that Li et al. [27] 
proves that the upper bound obtained by Corollary 3.2 is optimal. Moreover for some 
Gaussian anisotropic random fields, Xiao [37] also obtains a sample path regularity in 
the stronger L^-sense on whole the compact Kd- This follows from an extension of the 
Garsia-Rodemich-Rumsey continuity lemma Garsia et al. [16] or the minorization metric 
method of Kwapieh and Rosinski [22] . This would be interesting to study if these results 
still hold when considering a quasi-metric p (and not a metric) and if they can be applied 
to obtain the sample path regularity of S in the stronger L^-sense on whole the compact 
Kd, strengthening the assumption on the parameter s. 

3.2. Rate of almost sure uniform convergence 

This section is concerned with the rate of uniform convergence of the series (>S'Ar)^gfj 
defined by (4). Under Assumption 1, this series converges to S and, for any integer N, 
we consider the rest 


+ CO 

Rn{x) = S{x) - Sn{x) = ^ Wn{x)gn, XGKdCR^. 

n=N+l 


Then, conditionally on {Wn)n>i^ Rn{x) — Rn^v) is a sub-Gaussian random variable with 
parameter 


rN{x,y) 


f +00 '' 

^ \W^{x)-W^{y)t 

\n—N-\-l / 


1/2 


x,2/ G Kd. 


(9) 


Observe that Rq = S and that ro(x, y) = s(x, y). To obtain a rate of uniform convergence 
for the sequence the general assumption relies on a rate of convergence for the 

sequence (rwlTveN- 


Assumption 3. Let xq G Kd with Kd = C and let p be a quasi-metric 

on satisfying (7). Assume that there exist an almost sure event Lt', some random 
variables 7 > 0, /3 G (0,1], ry G R and a positive random sequence (&(iV))ArgN such that on 

n', 


WN €N,yx,y€B{xo,-f)r\Kd, rNix,y)<b{N)p{x,yf\log{p{x,y))\"’. (10) 

Note that fl', p and the random variables 7 ,/?,?? and h{N) may depend on a:o. Note 
also that since Assumption 3 implies Assumption 2, according to Theorem 3.1, almost 
surely, there exists 7 * G ( 0 , 7 ) such that i?w = S — Sm is continuous on B{xo,j*). The 
following theorem precises the modulus of continuity of Rn with respect to N and a rate 
of uniform convergence. 
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Theorem 3.3. Assume that Assumptions 1 and 3 are fulfilled. 

1. Then, almost surely, there exists 7 * € ( 0 , 7 ) and C € (0,+cx)) such that for 

\Rn{x) -RN{y)\ < Cb{N)y/log{N + 2)p{x, y)^\logp{x, 

for all N and all x,y G B{xo,"f*) fl Kd- 

2. Moreover, if almost surely, for all N gN, 

\RN{xo)\<b{N).y\og{N + 2), (11) 

then, almost surely, there exists 7 * G (0, 7 ) and C G (0, +c») such that 

\RN{x)\<Cb{N)./logiN + 2) 
for all N gN and all x G B{xo,j*) fl Kd- 

Proof. See Appendix B.2. □ 

An analogous of Corollary 3.2 holds for strengthening the previous local theorem to 
get uniform results on Kd or on B{xo,"f) n Kd when 7 is deterministic. 

Corollary 3.4. Assume that Assumptions 1 is fulfilled. 

1. Assume that Assumption 3 holds for any xq G Kd with the same almost sure event 
Tl', the same random variables (3 and rj, the same sequence ( 6 (A^))ArgN o,nd the same 
quasi-metric p. Then assertion 1 of Theorem 3.3 holds replacing B(xo,j*) Cl Kd by 
all the set Kd- If moreover, equation (11) is fulfilled for some xq G Kd, then 

7 7 

sup sup - , ^ < +00 almost surely. 

NefixeKd b{N)y/\og{N + 2) 

2. Assume now that Assumption 3 holds with a deterministic 7 . Then Theorem 3.3 
holds replacing B{xo,j*)r]Kd by B(xo,^) D Kd. 


4. Shot noise series 

4.1. Preliminaries 

In this section, we consider the sequence of shot noise series defined by 

N 

VIV G N, Va G = [a, 6] c (0,2), S% (a) = ^ A„, 

n—1 

where for all n > 1, the random variable T„ is the nth arrival time of a Poisson process 
with intensity 1 and (A„)„>i is a sequence of i.i.d. symmetric random variables, which 
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is assumed independent of (T„)„>i. Let us first recall that S*j^{a) converges almost 
surely to S'* (a) an a-stable random variable as soon as (see, [35], e.g.). Under 

a strengthened assumption on the integrability of rate of pointwise almost sure 
convergence and rate of absolute convergence have also been given in Theorems 2.1 and 
2.2 of [ 11 ]. 

Since (^ra)„>i rnay not be a sequence of sub-Gaussian random variables, we cannot 
apply Section 3 to the sequence (S^)^gpj. However, due to symmetry of (^n)„>i, 




n>l 


with (gra)„>i a Rademacher sequence independent of and Section 3 allows 

to study 

N 

SN{a) = ^ Wn{a)gn with Wn{a) := 
n—1 

Moreover, in Theorems 3.1 and 3.3, Sat (resp., Rjq = S — Sn) can be replaced by S'ff 
(resp., = S* — S^), see the proof of next theorem for details. Then, assuming that 
Xn is sufficiently integrable, we obtain the uniform convergence of S'f; on a deterministic 
compact interval Ki = [a, 6 ] C (0,2) and a rate of uniform convergence. These results, 
stated in the following theorem, strengthen Theorem 2.1 of [11] which deals with the 
pointwise rate of convergence. 


Theorem 4.1. For any integer n > 1, let he the nth arrival time of a Poisson process 
with intensity 1. Let {Xn)n>i be o, sequence of i.i.d. symmetric random variables, which 
is assumed independent of (r„)„>i. Furthermore assume that E(|Xip^’) < +cxd for some 

p>0. 

1. Then, almost surely, for all (0, min(2,2p)) and for all aG (0,6], the sequence of 

partial sums converges uniformly on [a, 6]. 

2. Moreover, almost surely, for all b G (0, min(2, 2p)) and for all a G (0, b), for all p' > 0 
with 1/p' G (0,1/6 — 1/min(2p, 2)), 


sup sup nGp' 

iVGN Q:G[a,&] 


+ 00 




< + 00 . 


Proof. See Appendix C.l. 


□ 


4.2. Modulus of continuity and rate of convergence of shot noise 
series 

In this section, we focus on some shot noise series, which are particular examples of condi¬ 
tionally sub-Gaussian series. For this purpose, we assume that the following assumption 
is fulfilled. 




12 


H. Bierme and C. Lacaux 


Assumption 4- Let (T'ra)„>i, (14,)„>i and (ffn)„>i be independent sequences satisfying 
the following conditions. 

1. {gn)n>i is a sequence of independent complex-valued symmetric sub-Gaussian ran¬ 
dom variables with parameter s = 1. 

2. Tn is the nth arrival time of a Poisson process with intensity 1. 

3. (T4,)„>i is a sequence of i.i.d. complex-valued random fields defined on Kd+i C 
(0,2) X 

4. For any (a,u) G LCd+i, Vn(a,u) G 

For any integer n > 1, we consider the complex-valued random field Wn defined by 

Wn(a,u):=T-^/‘^Vn(a,u), (a,u) G Ka+i C (0,2) x (12) 

Since |y„(a,u)|^ G and a/2 G (0,1), according to Theorem 1.4.5 of [35], 

+ CO +00 

\Wn(a,u)\'^ = y^T“^/“|l/„(a,M)|^ < -l-oo almost surely. 

n—1 n—1 

Therefore, the independent sequences (Wn)n>i and (gn)n>i satisfy Assumption 1. Then, 
S and (5'Ar)ArgN are well-defined on K^+i C (0,2) x C by (5) and (4). Before we 
study, the modulus of continuity of S and the rate of convergence of (Sn)^^^, let us 
state some remarks. 

Remark 4-1- Assume that conditions 1-3 of Assumption 4 are fulfilled with ( 3 n)„>i 
a sequence of i.i.d. random variables. Then Remark 2.6 of [34] proves that condition 4 
is a necessary and sufficient condition for the almost sure convergence of (S'Ar(a, M))^gjij 
for each (a,u) G AT^+i. Note that by Ito-Nisio theorem (see, e.g.. Theorem 6.1 of [23]), 
it is also a necessary and sufficient condition for the convergence in distribution of the 
sequence {Sn(o:,u)) Then, condition 4 is not a strong assumption and is clearly 
essential to ensure that S(a,u) is well-defined. 

Remark 4-2. Assume that Assumption 4 is fulfilled with (gn)„>i a sequence of i.i.d. 
random variables. Then, it is well known that for each a G (0,2), S(a,-) is an a-stable 
symmetric random field, as field in variable u. In Section 5.1, we will focus on a-stable 
random fields defined through a stochastic integral and see that, up to a multiplicative 
constant, such a random field Xa has the same finite distributions as S(a, ■) for a suitable 
choice of (gn,Vn)n>i- The sample path regularity of S in its variable a is not needed to 
obtain an upper bound of the modulus of continuity of Xa- Nevertheless, this regularity 
is useful to deal with multistable random fields (see Section 5.2). 

The sequel of this section is devoted to simple criteria, based on some moments of 14,, 
which ensure that Assumption 3 (and then Assumption 2) is fulfilled. More precisely, the 
results given below help us to give simple conditions in order to get Assumption 3 and 
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(11) satisfied with h{N) = {N + 1) for some convenient p' > 0. Then, all the results 
of Section 3 hold. 


Theorem 4.2. Assume that Assumption 4 is fulfilled with Kd+i = [a, b] x 11^=1 ^j] ^ 

(0,2) X and let p be a quasi-metric on R'^^^ satisfying equation (7). Assume also that 
for some Xg € itld+i, there exist r € (0, +oo), /3 G (0,1], 77 G R and p G (6/2, + 00 ) such 
that E(|Vi(a;o)|^^) < +00 and 


E 


sup 

x,yGKd+i 

0<||a;-j/||<r 


\Viix)-V,iy)\ 
p{x,y)P\logp{x,y)\'^ _ 


2p 


< 00 . 


(13) 


Let us recall that S and Sn are defined by (5) and (4) with Wn given by (12). 

1 . Then, almost surely (S'Ar)^gp,j converges uniformly on Kd+i and its limit S belongs 
almost surely to ?^p(it'd+i,/3,max( 77 , 0 ) + 1 / 2 ). 

2. Moreover, when p' > 0 is such that 1/p' G (0,1/6 — 1/min(2p, 2)), almost surely 

sup sup |«S'(a;) — 5'Ar(a;)| <+ 00 . 

Nen xGKd+i 


Proof. See Appendix C.2. □ 

Example 4-1 ■ Assume that Vi is a fractional Brownian field on R^^ with Hurst parameter 
H. Then (13) is satisfied for all p > 0 with p{x,y) = ||x — 7 /||, (3 = H and 77 = 1/2 (see, 
e.g., Theorem 1.3 of [ 6 ]). 

Let us now present a method (similar to those used in [7, 8 , 20] to bound some condi¬ 
tional variance) to establish (13). 

Proposition 4.3. Let xg = (ao,Mo) G A'^+i with = [a,b] x ]/[^_ ^i[a„ 6 ,] C (0,2) X 
R'^. Let Vi be a complex-valued random field defined on K^+i. Assume that there exists 
a random field (fy( 6 ))^gjg with values in [ 0 ,oo) and such that 

(i) there exists p a quasi-metric on R'^+^ satisfying equation (7) such that almost 
surely, 

Vx, 7/G ATd-Hi, 1 ^ 1 ( 2 ;) -Vi{y)\ <g{p{x,y)); 

(ii) there exists hg G (0,1] such that almost surely, the function h i-G g{h) is monotonic 
on [ 0 , 60 ]; 

(iii) there exist p > 6/2 and some constants /3 G (0,1], 77 G R and C G (0, 00 ) such that 
for some e > 0 and for h> 0 small enough. 


I{h) := W.{g{hfP) < Ch^P^\\og 7i|2p('7-i/2p-+ 


( 14 ) 





14 


H. Bierme and C. Lacaux 


Then, equation (13) holds for r > 0 small enough. 

Proof. See Appendix C.2. □ 

Remark 4-3. If {Vn)n>i is a sequence of independent symmetric random variables, 
Theorem 4.2 still holds replacing SN{ct,u) (resp., S{a,u)) by 

N / +00 

5;^(a,u) = ^r-i/“K(a,u) resp., by u) = ^ u) 

n—1 \ n—1 

In particular, following Example 4.1, when Assumption 4 is fulfilled with Vi a fractional 
Brownian field on with Hurst parameter H, assumptions of Theorem 4.2 are fulfilled 
with pd+i the Euclidean distance on (3 = H and rj = 1/2, on any compact {d + 1)- 

dimensional interval ATd+i- Especially, this leads to an upper bound of the modulus 
of continuity of S* on any compact {d + l)-dimensional interval Kd+i- Then for any 
fixed oq S ( 0 , 2 ), we also obtain that the oo-stable random field {S*{aQ,u)).,^^g^d is in 
'Hpj{Kd,H, 1) for pd the Euclidean distance on R'^ and for any compact set Kd C R'^. 

4.3. LePage random series representation 

Representations in random series of inhnitely divisible laws have been studied in [24, 25] . 
Such representations have been successfully used to study sample path properties of some 
symmetric oo-stable random processes {d = 1) and fields (see, e.g., [7, 8 , 12, 20]). 

Let us be more precise on the assumptions on the LePage series under study. 

Assumption 5. Let (Tn)„>i and {gn)n>i as in Assumption 4- Let (Cn)„>i be a 
sequence of i.i.d. random variables with common law 

equivalent to a a-finite measure v on (M'^, S(R‘^)) (that is such that > 0 for v-almost 
every (,). This sequenee is independent from (5n,7ra)„>i- Moreover, we consider 

Vn{a,u) := fa{u,fn)m{^n)~^^°‘, 

where for any a G Ki C (0,2), fa'. Kd X R'^ is a deterministie function such that 
VuG KdCM.‘^, [ 1/0(14,^)|“i^(d^) <+ 00 . 


Under this assumption. Assumption 4 is fulfilled with Kd+i = Ki x Kd. Then, empha¬ 
sizing the dependence on the function m, 5m, at and Sm are well defined on Kd+i by (4) 
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and (5) with Wn given by (12). In particular, 


Sjn{a,u) = 


(a, u) € Xd+i := KixKdC ( 0 , 2 ) x 


(15) 


Under appropriate assumptions on fa and m, the previous sections state the uniform 
convergence of the series, give a rate of convergence and some results on regularity for 
Sm- Precise results on regularity of Sm may be obtained using the following proposition, 
which states that the finite distributions of Sm does not depend on the choice of the 
i^-density m. 

Proposition 4.4. Assume that Assumption 5 is fulfilled and let Sm be defined by (15). 
Let (fn)n>i be a sequence of i.i.d. random variables with common law = TO(^)^(d^) 

equivalent to v. Assume that the sequences {f,n)n>i, (5n)„>i and (T'n)„>i are indepen¬ 


dent. 


1. Then, Sm^= Sm, where means equality of finite distributions. In other words, 



2. Assume moreover that for v -almost every ^ the map (a, u) i—>■ /q(u, ^) is con¬ 

tinuous on the compact set K^+i C (0,2) x Let us consider p a quasi-metric on 
/3 G (0,1] and 77 G R. Then, Sm belongs almost surely in TLp{Kd+i,fi, rf) if and 
only if Sm does. 


Proof. See Appendix C.3. 


□ 


In particular, when studying the sample path properties of Sm, this result allows us 
to replace m by an other function rh so that the regularity of Sm may be deduced 
from the regularity of Sm- For example, replacing m by depending on a:o this may 
lead to a more precise bound for the modulus of continuity of Sm around xq (see, e.g., 
Example 5.3). 

5. Applications 

5.1. CK-stable isotropic random fields 


Let us fix a = oq S (0,2) and assume that Assumption 5 is fulfilled with pn some isotropic 
complex random variables. Then, the proof of Proposition 4.4 (see Section C.3) allows to 
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compute the characteristic function of the isotropic oo-stable random field 5'm(Q!or) = 
iSmiao,u))ueK^’ which leads to 

Sm{ao, •) '= ( [ fao{u,OMc,MO) 

with Mao ^ complex isotropic ao“Stable random measure on with control measure ly 
and 

/ 1 \ / /* + 00 Q* (fi\ \ 

<=E(|JR(gi)r“)'/“”(^—|cos(0 )rM0j ^ ■ (^6) 

When z/ is a finite measure (resp., the Lebesgue measure), this stochastic integral repre¬ 
sentation of Sm{cio,-) has been provided in [29, 35] (resp., [7, 20]). 

Let us note that assumptions of Theorem 4.2 and Proposition 4.3 can be stated in term 
of the deterministic kernel fao to obtain an upper bound of the modulus of continuity of 
Sm- In general, well-choosing and applying Proposition 4.4, we obtain a more precise 
upper bound of the modulus of continuity of Sm{aQ, ■) around mq, which also holds for a 
modification of the random held 

Xao=([ faoiu^mcaido) ■ (17) 

VJR'^ / ttGifd 

To illustrate how the previous sections can be applied to study the Held Xao, which 
is dehned through a stochastic integral and not a series, let us focus on the case of 
harmonizable stable random Helds. More precisely, we consider 

/„„(u,e) = (e‘<“’«>-!)./>«„ (0, V(u,e)GM''xR^ (18) 

with ipao : —)• C a Borelian function such that 

/ min(l,lieir“)|z/>ao(e)rV(de)<+oo. 

JR'^ 

Note that, since this assumption does not depend on u, the random Held Xao may be 
dehned on the whole space For the sake of simplicity, in the sequel, we consider 
the case where v is the Lebesgue measure and Hrst focus on a random Held Xao which 
behaves as operator scaling random Helds studied in ]9]. 

Proposition 5.1. Let Oq G (0,2) and let Xao defined by (17) with v the Lebesgue 
measure on R"^. Let E be a real matrix of size d x d whose eigenvalues have positive 
real parts. Let te and te* be functions as introduced in Example 2.1 and let us set 
q{E) = trace(£’) and ai = min;,ggp(£;) 5R(A) with Sp{E) the spectrum of E, that is, the 
set of the eigenvalues of E. Assume that there exist some finite positive constants c.^, A 
and (3 G (0,ai) such that 

IV'ao(C)l < for almost every U^]] > A. (19) 
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Then, there exists a modification of Xag such that almost surely, for any £ > 0, for 
any non-empty compact set Kd C 


sup 


\X*{u)-X*iv)\ 


te{u — u)^[log (1 + te{u — v) ^ ^ ^ ° 

U^V 


< +CX). 


Remark 5.1. The quasi-metric {x,y) te{x — y) may not fulfill equation (7) since 
the eigenvalues of E may not be greater than 1. Nevertheless, the quasi-metric {x,y) i—> 
TE/aiix — y) does and the conclusion with te in the previous proposition then follows 
from the comparison 


e CiTEifT^ < ^E/ai iO < C2TE{fT^ 

with Ci,C 2 two finite positive constants. 

Proof. See Appendix D.l. □ 

An upper bound for the modulus of continuity of such harmonizable random fields 
is also obtained in [38]. This upper bound is given in term of the Euclidean norm and 
then does not take into account the anisotropic behavior of Xag. Even when te is the 
Euclidean norm, our result is a little more precise than the one of [38]. The difference is 
only in the power of the logarithmic term. 

Let us now give some examples. We keep the notation of the previous proposition and 
the eigenvalues of the matrix E have always positive real parts. 

Example 5.1 (Operator scaling random fields [9]). Let ^ [0,oo) be an E*- 

homogeneous function, which means that 

Vce (0,-i-oo),V^e]R‘', f) = 

where c^ = exp (A* log c). Let us assume that ip is a continuous function such that 
tpif) ^0 for ^ 7 ^ 0. Then we consider the function tpag : —>■ [0, -foo] defined by 

The random field Xag, associated with tpag by (17) and (18), is well-defined and is 
stochastically continuous if and only if 77 S (0,ai). Then, let us now fix 77 g (0,ai). 
Since "ipag is 7f‘-homogeneous, one easily checks that there exists c.^ € (0, -foo) such that 

Ve G V’ao (?) < (e)-^-^(^)/“» . 

Then, the assumptions of Proposition 5.1 are fulfilled with /3 = 77. The corresponding 
conclusion was stated in Theorem 5.1 of [7] when 77 = 1 and ai > 1, which is enough to 
cover the general case using Remark 2.1 of [7]. 
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Example 5.2 (Anisotropic Riesz Bessel cx-stable random fields). Let us con¬ 
sider 


V'ao (■?) 


1 


with two real numbers /3i and /32- Assuming that 


eeR''\{ 0 } 


QjE) 

2 


< Pl+ 132 


and 



OfOOi 

2 


the random field Xa^ is well-defined by (17). When is the Euclidean norm, this 
random field has been introduced in [38] to generalize the Gaussian fractional Riesz- 
Bessel motion [2]. 

We distinguish two cases. If /3i -1-/32 < Proposition 5.1 can be applied with 

P _ ‘^{Pi+P 2 )^-q{E) ^ Proposition 5.1 can be applied for any (3 S (0,ai). 


Random fields defined by (17) have stationary increments so that their regularity on 
Kd does not depend on the compact set Kd- To avoid this feature, one can consider 
non-stationary generalizations by substituting ipao by a function that also depends on 
u € Kd ■ More precisely, we can consider 


Xc.o=( f (eb“-«)-l)^„„(w,e)M„„(do') (20) 

\dR<i / u^Kd 

with Mag a complex isotropic ao-stable random measure with Lebesgue control measure 
and tpao 8. Borelian function such that, for all u S Kd, 

f |eb“’«> - l|“”|^/>„„(M,^)|“”dC < +00. 

Under some conditions on 'ipag, when considering the local behavior of Xag around a 
point uq one can conveniently choose a Lebesgue density to obtain an upper bound 
of the modulus of continuity of the shot noise series {ao, •) given by (15) with 

- l)V'ao(W:0- 

For the sake of conciseness, let us illustrate this with multi-operator random fields, which 
have already been studied in [8] . 


Example 5.3 (Multi-operator scaling cx-stable random fields). In [8], we consider 
E a function defined on with values in the set of real matrix of size dx d whose eigen¬ 
values have real parts greater than 1 and tjj x ^ [0,-|-oo) a continuous function 
such that for any u € R^^, i){u, ■) is homogeneous with respect to E{uy, that is, 

'0(m,c'®^“^ ^) = c^('u,^), V^gR'^,Vc>0. 
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Under convenient regularity assumptions on if and E, the ao-stable random field Xa^ is 
well-defined by (20) setting 

ifaa ^)-l-9(-E(u))/ao q[E{u)) = traCe(i?(M)). 

Let Kd = 0^=1 C and uq £ Kd- Let us set Kd+i = {ao} x Kd and consider the 

quasi-metric p defined on R'^+^ by 


p((a,u), (a',u)) = \a- a '\-t T£;(„„)(m - v) 


for all (a,u), {a',v) £ M x R'^, which clearly satisfies equation (7). Then, under assump¬ 
tions of [8] , there exists a Lebesgue density m„(, > 0 a.e. such that Assumption 2 holds 
for on Kd+i with 77 = 0 and all fd £ (0,1), adapting similar arguments as in Proposi¬ 
tion 5.1 (see Lemma 4.7 of [8]). Therefore, following a part of the proof of Proposition 5.1, 
there exists a modification X*^ of X^f, such that almost surely. 


lim sup 

^ 1 ^ u,veiB{uQ,r)nKfi 
u^v 


- Xl^{v)\ 

Te{uo){u - 


< -|-oo 


for any e £ (0,1). This is Theorem 4.6 of [8]. 

For the sake of conciseness, we do not develop other examples. Nevertheless, let us 
mention that our results can also be applied to harmonizable fractional a-stable sheets 
or even to operator stable sheets. In particular, this improves the result stated in [30] for 
fractional a-stable sheets. Note that we can also deal with real symmetric measure Wa- 

5.2. Multistable random fields 

Multistable random fields have first been introduced in [13] and then studied in [14]. 
Each marginal X{u) of such a random field is a stable random variable but its stability 
index is allowed to depend on the position u. 

Generalizing the class of multistable random fields introduced in [26], we consider a 
multistable random field defined by a LePage series. More precisely, under Assumption 5, 
we consider 


+ 00 

= ueKd, ( 21 ) 

n— 1 

where a : Kd —>■ (0,2) is a function. Then since Sm{u) = Sm{o:{u), u) with Sm defined by 
(15), we deduce from Section 4 an upper bound for the modulus of continuity of S. In 
particular, assuming that a is smooth enough, we obtain the following theorem. 
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Proposition 5.2. Let Kd = 11^=1 ^ Let us choose uq S it'd- Let p be a quasi¬ 

metric on satisfying equation (7) and let a:Kd^ (0)2) belongs to ?ip(it'd, 1,0). Let 
us set 

a = min a, & = maxa and isTi = [a, &] C (0,2) 

Kd Kd 

and consider the quasi-metric p defined on R x R"^ by 

p{{a, u), (o', n)) = la — a'l + p{u,v). 

Assume that Assumption 5 is fulfilled and that equation (13) holds on Kd+i = [a,b] x it'd 
for some p > 6/2, /3 £ (0,1] and 77 £ R. Assume also that 

E(|Pi(a(uo),Mo)|^^) = / <+00. 

dR<i 

Let Sm,N 6e defined by ()) with Wnicx^u) = Tn faiu,f,n)'m{f,n) and let Sm,N{u) = 
Sm,Nia{u),u). 

1. Then, almost surely, {Sm,N)N£n converges uniformly on Kd to Sm and almost 
surely the limit Sm belongs to TLpiKd, I3,m.a,yi{p,0) + 1/2). 

2. Moreover, for all p' > 0 such that 1/p' £ (0,1/6 — 1/min(2p, 2)), 

sup sup |,Sm(M) - Sm,N{u)\ < +00. 

Nen ui^Kd 

Proof. See Appendix D.2. □ 

Remark 5.2. Let us recall that Sm € iip(ifd,/?, max ( 77 ,0) + 1/2) if and only if Sm € 
'Hp{Kd,l3, max ( 77 ,0) +1/2), with rh an other z 7 -density equivalent to v, by Proposition 4.4. 

To illustrate the previous proposition, we only focus on multistable random fields 
obtained replacing in a LePage series representation of an harmonizable operator scaling 
stable random field the index a by a function. Many other examples can be given, such 
as multistable anisotropic Riesz-Bessel random fields or the class of linear multistable 
random fields defined in [14]. 

Corollary 5.3 (Multistable versions of harmonizable operator scaling random 
fields). Let E be a real matrix of size dxd such that minAgSpB 3i(A) > 1. Let us consider 
Pe and te as defined in Example 2.1. Let us also consider —>■ [0,oo) a continuous, 

E'^-homogeneous function such that 7 ^ 0 for Then we set 

/„(n,e) = (e‘<“’«>-l)7/(e)-i-«(^)/“ 

with q{E) =trace(i?). Let m be a Lebesgue density a.e. positive on M.'^, {f,n,Tn,9n)n>i 
as in Assumption 5 with v the Lebesgue measure and consider a function aiR"^ —>■ (0,2). 
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Therefore, the multistable random field Sm is well-defined by (21) on the whole space 
Moreover if a & ,!,{)), then for anyuo^W^ and£>0, there exists r € {0,1] such 

that almost surely 


sup - 

u,veB(uo,r) Te(u 

U^V 


v)\\OgTE{u — + 


< +CX). 


Proof. See Appendix D.2. □ 

Remark 5.3. In particular, when E = Id, te is the Euclidean norm and we obtain an 
upper bound of the modulus of continuity of multistable versions of fractional harmoniz- 
able stable fields. 


Appendix A: Proof of Proposition 2.1 

The proof of Proposition 2.1 is based on the following lemma. 

Lemma A.l. If Z is a complex-valued sub-Gaussian random variable with parameter 
s £ (0, +oo), then for all t £ (0, +oo), P(|^| > t) < 4e“‘ ^. 

Proof. Let t £ (0,+oo). Since Z is sub-Gaussian with parameter s, and are 
real-valued sub-Gaussian random variables with parameter s. Then applying Proposi¬ 
tion 4 of [17], 

P(|Z|>t)<p(^|K(Z)|>0 +p(^|cj(Z)|>0 <4exp(^-^), 
which concludes the proof. □ 


Let us now prove Proposition 2.1. 

Proof of Proposition 2.1. Let t £ (0,-|-oo). Since Proposition 2.1 is straightforward 
if a = 0, we assume that a ^0. Since the sequence (gn)„>i is symmetric, by the Levy 
inequalities (see Proposition 2.3 in [23]), for any M £ N \ {0}, 


sup 

P 

\l<P<M 

n—1 


>t||a||^2 <2P 


M 


n—1 


>t a 




We now prove that is sub-Gaussian. By independence of the random variables 

Pn and since each pn is sub-Gaussian with parameter s = 1, 


M 


M 


Vz£C, 


E(e 






,a?(za„g„) 


)<n 




( 22 ) 
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with SM = {Yjn=i ^ ll«ll^ 2 . Hence, for any M € N \ {0}, Y.n=i<^ngn is sub- 

Gaussian with parameter sm- Since a ^ 0, for M large enough, sm ^ 0 and then applying 
Lemma A.l, 


Vt >0, 


sup 

p 

^ ^ ^nQn 

y i<p<M 

n—1 



< 8 exp 




< 8 e"*"/®. 


Assertion 1 follows letting M —)■ +oo. 

Let us now prove assertion 2. If there exists A^ G N \ {0}, such that 

Vn >N, an = 0, 

then, assertion 2 is fulhlled since o-ng-n = is a sub-Gaussian random 

iV 2 ^ 

variable with parameter sjv = (Sn=i l®"l ) = 11^11^2■ Therefore to prove assertion 2, 

we now assume that 

yN GN\{0},Bn>N, an^O, 

so that \o.n\^ 7^0 for any integer A^ > 1. Then, applying assertion 1 replacing an 

bya„l n>Af, we have 


Ve>0,VAf gN\{0}, 


sup 

P 

E °‘^9n 

\P>N 

n^N 





Since ||a ||^2 = |on|^ < +oo, this implies that (X]^=i o.ngn)j^ is a Gauchy sequence in 

probability. Then, by Lemma 3.6 in [18], the series J2n=i^ngn converges in probability. 
By It6“Nisio theorem (see, [23], e.g.), this series also converges almost surely, since the 
random variables n > 1, are independent. Moreover, since sup^>]^ = Uoj]^^ < +oo, 

equation (22) implies the uniform integrability of the sequence for 

any z gC. Then, letting M -> +oo in (22), we obtain that an5n is sub-Gaussian 

with parameter l|aj|^ 2 . Moreover, we conclude the proof noting that 


Vt>0, 


+ 00 

E 

n—1 




>f||a|| 


t2 I < 


sup 
. P>1 


E 

n—1 


O^nQn 


>t|]c 


lh2 ^ 


<8e-‘'/s. 


□ 


Appendix B: Main results on conditionally 
sub-Gaussian series 

B.l. Local modulus of continuity 

This section is devoted to the proofs of the results stated in Section 3.1. 
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Proof of Theorem 3.1. Let us recall that xq € Kd = Y[j=i['^jjbj] C We assume, 
without loss of generality, that 


VI < j < d, Oj < bj. 

Actually, if some aj = bj, we may identify limits S as random fields 

defined on Kd' C for d' < d. Note that if aj = bj for all 1 < j < d, there is nothing to 
prove. 

We also assume that 7 ( 0 ;) € (0,1), which is not restrictive and allows us to apply 
equation (7) as soon as ||x — y\\ < 7 ( 0 ;) (with C 2 ,i and 02,2 which do not depend on 7 ). 

First step. We first introduce a convenient sequence of countable sets included 

on dyadics, which is linked to the quasi-metric p. It allows to follow some arguments of 
the proof of the Kolmogorov’s lemma to obtain an upper bound for the modulus of 
continuity of S. 

Let us hrst introduce some notation. For any fc S N \ {0} and j G Z^, we set 

Xk,j = -^, 'Dk = {xk,j-.j and = min{n e N \ { 0 } : 02,2^—< 2 "'"} 

with C 2,2 the constant given by equation (7). Then, choosing 02,2 large enough (which 
is not restrictive), one checks that {vk)k>i i® increasing sequence. In particular, the 

sequence (77i,^)^>^ is increasing and = UiS = US ■ Moreover, V f] Kd is 
dense in Kd since Oj < bj for any 1 < j < d. Then, as done in Step I of the proof of 
Theorem 5.1 of [7], one also checks that for k large enough, 'D,y^. n Kd is a 2“^ net of Kd 
for p, which means that for any x G Kd, there exists j G Z'^ such that p{x,Xy,^j) < 2~^, 
with J =31^'"'" S Kd- 

Second step. This step is inspired from Step 2 of [7, 8 ]. The main difference is that we 
use Proposition 2.1 to obtain a uniform control in N. 

For k gN \ {0} and {i,j) G Z'^, we consider 

Elj = Iw: sup|S'Ar(a;^,,,) - 5Ar(a;i.,7)| > 

'■ Nm '* 

with, following [19], 

ipit) = y 8 Ad log i, t > 0, (23) 

for A > 0 conveniently chosen later. We choose S G (0,1) and set for fc G N \ {0}, 

4 = 

(24) 

4 = {(*, j) G (Z'^ n2‘''“A:d)^ < 4} and Ek= |J E^j. 

Since (/j is a decreasing function and s > 0, for any k gN \ {0} and for any (f, j) G Ik 

P(F;'] ) < P ( sup I S'at (xj.^ ,i) - S'at (xj.^ J) I > s (x^, J-) (^(4)) ■ 

'-Nen ' 
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Since {gn)n>i is a sequence of symmetric independent sub-Gaussian random variables 
with parameter s = 1, conditioning to (W^n)„>i and applying assertion 1 of Proposi¬ 
tion 2.1, one has 

Vfc G N \ {0}, V(i, j) G Ik, 

by definition of s, Sn, <P and 6k- Moreover, since Kd C is a compact set, using equation 
(8) and the definition of Vk, one easily proves that there exists a finite positive constant 

Cl G (0, -foo) such that for any fc G N \ {0}, card/fc < Hence, 

+ 00 +CXD +00 

HEk) <Y ^ Cl ^ ^-(A(l-S)-2/H+{l-S)/H)kdlos2 ^ 

k=l k=l(i,j)elk k=l 

choosing A> ^ — = and 6 small enough. Then, setting 


C + OO +CX) 

u 

k=lt=k 

with the almost sure event introduced by Assumption 2, the Borel-Cantelli lemma 
leads to P(H") = 1. Moreover, by Assumption 2, for any w G O" there exists k*{uj) such 
that for every k > k*{uj) and for all x,y € with x,y G B{xq,^{lo)) n Kd and p{x, y) < 
4 = 


sup |5'7v(a;) -5'Ar(y)| < Cp(a;,?/)^|log(p(a;,?/))|’'+^^^. (25) 

iVeN 

Third step. In this step, we prove that (25) holds, up to a multiplicative constant, 
for any x,y GV closed enough to xq. This step is adapted from Step 4 of the proof of 
Theorem 5.1 in [7], taking care that (25) only holds for some x,y G n Kd randomly 
closed enough of xq. Let us mention that this step has been omitted in the proof of the 
main result of [8] but is not trivial. We then decide to provide a proof here for the sake 
of completeness and clearness. 

Let us now hx w G H" and denote by k > 1 the constant appearing in the quasi-triangle 
inequality satisfied by p. We also consider the function F defined on (0, -boo) by 

F{h) :=/i'^|log(/i)r+^/l 

Observe that A is a random function since /? and p are random variables. Then, we 
choose ko = fco(w) G N such that the three following assertions are fulfilled: 

(a) F is increasing on (0,(5fe(,], where 6k is given by (24), 

(b) for all k > ko{Lu), 0 is a 2“^ net of Kd for p, 

(c) 2'=»4 o+i>3k2 . 
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Even if it means to choose k*{uj) larger, we can assume that k*(uj) > fcg and that 

/ r \ 1/H 

where H_ and C 2,2 are defined in equation (7). 

Let us now consider x,y GT> Ci Kd n i?(a;o, 7 *(a;)) such that x ^y. Let us first note 
that x,y & i?(a;o, 7 (w)). Moreover, since \\x — y\\ < 2'^*{uj) < 7 ( 0 ;) < 1, the upper bound 
of equation (7) leads to 


3K'^p{x,y) < 3 k‘^C2,2\\x - y\\— < 4 *(<^) 

by definition of 7 *(a;). Then, there exists a unique k > k*{u}) such that 

Sk-i-i <3K'^p{x,y) <6k. (27) 

Furthermore, since x,y GV Cl Kd, there exists n > fc + 1 such that x,y & fl Kd and 
for j = fc,..., n — 1 , there exist x^^'> S n Kd and y^^^ S T),,. n Kd such that 

p{x,x^^^) <2~^ and p{y,y^^^)<‘2~G (28) 

Let us now fix G N and focus on Sn{x) — SN{y)- Then, setting = x and = y, 

n — 1 

Sn{x) - SN{y) = ^ - Sn{x^^^)) 

j—k 

(29) 

n—1 

j=k 

The following lemma, whose proof is given below for the sake of clearness, allows to 
apply (25) for each term of the right-hand side of the last inequality. 

Lemma B.l. Choosing k*{w) large enough, the sequences and {y^^'^)k<j<n 

satisfy the three following assertions. 

{a.) x^^\y^^1 € B{xo,'y{ui)) for any j = k,...,n, 

(b) for any j = k,...,n-l, max(p(x(^+^\x('^)),p(y(-^+i),y(^))) < 4 + 1 , 

(c) p(x('=\y('=)) < 4- 

Therefore, even if it means to choose k*{ui) larger, applying this lemma and equations 
(25) and (29), we obtain 


/ n—1 

\Sn{x) - 5w(2/)| < C F(p(x('=), yW)) + 2 ^ ^(^.+ 1 ) 


j=k 
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since F is increasing on (0,(5^^] and since j >kQ. This implies, by definition of F that 

\Sn{x) - SN{y)\< C(F(p(xW,y('=))) + 2(7F(4+i)), 


C(w) = 2'^5f"\j + i)“a-('7(-)+i/2.o) < 
j=o 

since /3 > 0 and Sj = with ^ < 1. Then, since F is increasing on (0,4)j by 

assertion 3 of Lemma B.l and equation (27), we get 

|5iv(a:) - SNiy)\< C{1 + 2C)F{iK^p{x,y)), 

for every N and x,y GVn B(xo,"/'*(uj)) fl Kd- Therefore, by continuity of p and each 
SM and by density of 2? H Kd in Kd 

|5w(x) - SN{y)\< C{1 + 2C)F{3^p{x,y)), (30) 

for every N €N and x,y € B{xo,j*{u})) fl Kd- 

Fourth step: Uniform convergence of Sn- Let us now consider 

17 = [^ /S'Ar(w) = S'(m)| n n". 

Observe that P(f2) = 1. Let us now fix w G (7. Hence, by equation (30), the sequence 
(5'Ar(-)(<^))AreN> which converges pointwise on 27ni?(a;o, 7 *(a;)) OiGd is uniformly equicon- 
tinuous on i?(xo, 7 *(a;)). Since V D B{xo,'y*{uj)) CiKd is dense in H(a;o, 7 *(a;)) HKd, by 
Theorem 1.26 and adapting Theorem 1.27 in [33], converges uniformly on 

B{xo,^*{u:)) CiKd- Therefore, its limit S is continuous on B{xo,"t*{uj)) UKd- Moreover, 
letting N +oo in (30) (which holds since w € H), we get 

\S{x) - 5(2/)! < C(1 + 2C)F{3K^p{x,y)), (31) 

for every x,y G B{xo,j*{uj)) 0 Kd, which concludes the proof. □ 


Let us now prove Lemma B.l. 


Proof of Lemma B.l. Let us first observe that = x G B{xq,^{w)) UKd and = 
y G B(xq,^(u})) n Kd- Let us now fix j G {k ,..., n — 1}. The lower bound of equation (7) 
leads to 


- xjj + jjx - CColl < 


pO 


.U) 


,1/H 


^ 2,1 


— a^oll < 


+ ||a;-a;o||. 
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Since x G i?(xo, 7 *(w)) with 7 * satisfying equation (26) and since p{x^^\x) <2 ^ with 
j>k>k*iuj), 


— a;o|| < 


2 ”fc* / H 

^ 2.1 



Then, choosing k*{w) large enough, x^^'> € B(xo,j(uj)) for j = k,...,n— 1. The same 
holds for . Assertion 1 is then proved. 

Let us now observe that since j > fcp and since k > 1 , 


2J'5j+i >2'=“4o+i >3k2>3«; 


(32) 


by definition of fco (see the third step of the proof of Theorem 3.1). Then, using the 
quasi-triangle inequality fulfilled by p and (28), we obtain that 

p{x^^'^^\x^^^) < 3 k 2 “^'^+^^ < < 4+1 • 

Since the same holds for p{y^^^^\y^^')), assertion 2 is fulfilled. Moreover, applying twice 
the quasi-triangle inequality fulfilled by p and equations (27), (28) and (32) (with j = k), 
we obtain 

+p(a;,y)) < 3K^p{x,y) < 4, 

which is assertion 3. □ 


Let us now focus on Corollary 3.2. Its proof is based on the following technical lemma. 

Lemma B.2. Let = Ilfci 4] C be a compact d-dimensional interval, j5 G (0,1), 
7 G M and p be a quasi-metric on satisfying equation (7). Let {fn)n^f^ be sequence of 

o 

functions defined on Kd and let {B{xi,ri))^^^^^ be a finite covering of Kd by open balls 
with Xi G Kd and ri > 0. Assume that for each 1 <i <p, there exists a finite positive 
constant Ci such that 

\/x,y&B{xi,rf)nKd, sup|/„(a;) - fn{y)\ <Cip{x,yf\fog{l+p{x,y)~^y\'. 

riGN 

Then there exists a finite positive constant C such that 

Vx,yeKd, snp\fnix)-fn{y)\<Cp{x,yf[log{l-\-p{x,y)~^)]^. (33) 

raGN 


Proof. By the Lebesgue’s number lemma, there exists r > 0 such that 


yx e Kd,31 <i<p, B{x,r) cB{xi,ri). 
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Let us first note that since the map Fp : {u,v) !->■ z;)[log (1 + p{u,v) is positive 

and continuous on the compact set K = {{u,v) € K x iL/||u — z;|| > r}, 


m := inf Fp G ( 0 , +oo). 


K 


Then distinguishing the cases ||a: — y\\<r and jja; — y\\ > r, one easily sees that 



where M = sup^gpjsup^, \fn{x) — fniy)\- It remains to prove that M < +oo. Note 
that 


sup |/„(x) - fn{y)\<c max Cj, 

x,yGKfi 

l\x-y\\<r 


l<i<p 


where c = sup^g^^^^ Fp < +oo by continuity of Fp on the compact set x Kd- Then 
since Kd is a compact convex set, using a chaining argument, one easily obtains that 


M < + 00 , which concludes the proof. 


□ 


Proof of Corollary 3.2. We only prove assertion 1. Actually, assertion 2 is proved 
using the same arguments but replacing Kd by B{xo,'y) Kd- 

Assume that for any xq € Kd, Assumption 2 holds with D,', rj and the quasi-metric 
p independent of xq. Following the proof of Theorem 3.1 and keeping its notation, let us 
quote that 7 * and Q do not depend on xq- Let us now fix a; G 17. From the third step of 
the proof of Theorem 3.1 and Lemma B.2, we deduce that equation (30) still holds for 
any x,y G Kd- This allows to replace B{xo,'y*{uj)) by Kd in the fourth step of the proof 
of Theorem 3.1, which leads to assertion 1. □ 

B.2. Rate of almost sure uniform convergence 

Proof of Theorem 3.3. Let us first observe that Theorem 3.1 holds. Then, for almost oj, 
even if it means to choose 7 smaller, the sequence of continuous functions (<S'Ar(-)(a;))^gpj 
converges uniformly on i?(a;o, 7 ( 0 ;)) r\Kd, which implies that each Rn{-){uj) is continuous 
on B{xo,"f{uj)) nKd- As in the proof of Theorem 3.1, we assume without loss of generality 
that Kd = n^=iwith aj < bj. 

Proof of assertion 1. Since it is quite similar to the proof of equation (30), we only 
sketch it. 

For fcGN\{0}, iVeN and (z,j)eZ^, we consider 


Kf = - i?Ar(x^^ 7 )| > \/log(7V + 2)rN (x^^)V 3 (p(xj.,,i, x^^ 7 ))} 
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with rjv defined by (9), ip by (23) and by Step 1 of the proof of Theorem 3.1. 

Then, we proceed as in Step 1 of the proof of Theorem 3.1 replacing the set Ek by 


+ 00 




u u 


E: 


k,N 


with Ik and 6k defined by (24), and applying assertion 2 of Proposition 2.1 instead of 
assertion 1. Then, choosing the constant A, which appears in the definition of ip, and 
5 G (0,1) such that 

2 \ 3 

A{l — 6) — — -\ — :=->Q and A(1 — <5)log2 > 1 

we obtain that 

+ 00 + 00 +00 

'^f’iE'k) < C2 ^ = C2 X! < +00 

k=l N=2 N=2 

with C2 a finite positive constant. Then, by Borel-Cantelli lemma, the definition of ip and 
Assumption 3, almost surely there exists an integer k*(uj) such that for every k > k*(uj), 
for all G N, and for all x,y € with x,y € B{xq,^{u))) fl Kd and p{x,y) < 6k = 
2^ — {l—S)k 


\RNix) - RN{y)\ < Cb{N)^J\og{N + 2)p{x,yf\\og{p{x, j/)). 

In addition, replacing in Step 2 of the proof of Theorem 3.1, Sn by Rn (which still be, 
for almost all w, continuous on B{xq,^{uj)) ClKd), we obtain that for almost all w, there 
exists 7 * G ( 0 , 7 ), such that 

|i?jv(a:)-i? 7 v(?/)| <Cb{N)^/\og(I^^p{x,y)^\log{p{x,y))\'^^^^'^ (34) 

for every iV G N and x,y € B{xo,J*{uj)) n Kd- This establishes assertion 1. 

Proof of assertion 2. This assertion follows from equations (34) and (11), the conti¬ 
nuity of p on the compact set B{xo,j{u!)) DKd and 

|i?Ar(x)| < \Rn{x) - i?Ar(a;o)| -f |i?Ar(a;o)|. 

The proof of Theorem 3.3 is then complete. □ 

Proof of Corollary 3.4. We only prove assertion 1. Actually, assertion 2 is proved 
using the same arguments but replacing Kd by B{xq,^) H Kd- 

Let us assume that Assumption 3 holds with IV, (3, rj and the quasi-metric p inde¬ 
pendent of Xq. Note first that the almost sure event 11 under which (34) holds does not 
depend on xq. Then applying Lemma B.2 to /„ = i?„/ {b{n)y^log{n 2)), we obtain that 
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equation (34) still holds for x,y G Kd- If moreover for some xq, equation ( 11 ) is fulfilled, 
then following the proof of assertion 2 of Theorem 3.3, we also have: there exists C a 
finite positive random variable such that for all iV S N, 

sup \RN{x)\<Cb{N)^/\og(WT^, 
xGifd 

which concludes the proof. □ 


Appendix C: Shot noise series 

C.l. Proof of Theorem 4.1 

Let ( 5 n)„>i be a Rademacher sequence, that is, a sequence of i.i.d. random variables 
with symmetric Bernoulli distribution. This Rademacher sequence is assumed to be in¬ 
dependent of {Tn, Then, by independence and also by symmetry of the sequence 

(^n)„>i, {Xngn)n>i has the Same distribution as (^n)„>i and is independent of the 
sequence (r„)„>^. 

Let us now set 


N 

Wnia) =T~~^/°‘Xn and SN{a) ='^Wn{a)gn, 

n—1 

SO that {S^{a)^N > 1} has the same finite distribution as {SN{a),N > 1}. Moreover, 
since 

N N 

^|IT„(a)|" = ^r-2/“|X„p 

n—1 n—1 

with Xn G (with p > 0), Assumption 1 is fulhlled on any Ki = [a, b\ C (0, min(2, 2p)) 
(see, e.g., [35]). 

Let us now hx a, 6 G (0, min(2, 2p)) such that a <b, a' G (0, a) and b' G ( 6 , min(2, 2p)). 


Proof of assertion 1. By the Mean Value Inequality, we get that for any a^a' G [a,b] 
and n > 1 , 


< cjo - o'j max(T-i/'>',T-i/a') ( 35 ) 

with c a finite positive constant. It follows that, almost surely, for all a, a' G [a, 6 ], 



with C = c{Y,nTi Tn \Xn\’^ + \Xn\^Y^‘^ < +00 since |V„p e LP with 2p > 

b' > a' and a', b' G (0,2). Therefore, the assumptions of assertion 1 of Corollary 3.2 hold. 
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Let us now remark that for all a, a' G [a, 5], 

{(^■^(a) - s(a, a')); ^ > 1} '= s(a, a')); N > 1}. 


This allows us to replace Sn by 5'^ in the second step of the proof of Theorem 3.1. Then, 


the third and the fourth step of this proof still hold replacing Sat by S’f; and the limit S 


by the limit S* since each S*j^ is continuous (as Sn is) and since converges pointwise 
to S*. This allows us to also replace {Sn,S) by {S'ff,S*) in the proof of assertion 1 of 
Corollary 3.2. It follows that almost surely, (<S')^)jvgN converges uniformly on [a,h\ to S*. 
Since this holds for any 0 < a < b < min(2, 2p), assertion 1 of Theorem 4.1 is established. 

Proof of assertion 2. Since almost surely the sequence of continuous random fields 
converges uniformly on [a,&], for all G N the rest defined by 


RUc^)--= E 


n—N+1 


is also continuous on [a,b]. Remark also that we have, for all a,a' G [a,b] and iV G N, 


{R*N{a) - R%{a'),rNia,a')) = (RAr(a) - i?Ar(a'):^Af(«!«')) 
where i?Ar(a) = I]nrAr+i Tn^^°‘Xngn = S{a) - S'Ar(a) and 



As done for Sn, the previous lines allow to replace Rn by in the proof of Theorem 3.3. 
Moreover, by equation (35), almost surely, for all G N, and a,a' G [a,b], 


1/2 



rN{a,a')<c\a-a'\( E E 


(37) 






Let us now fix p' > 0 such that 1/p' G (0,1/6— 1/min(2p,2)). Choosing if necessary 
b' > b smaller, we assume without loss of generality that 1/p' G (0,1/6' — 1/min(2p, 2)) C 
(0,1/a' — l/min(2p,2)). Then, by Theorem 2.2 in [11], almost surely, for all a,a' G [a,6], 



since X^ G with p > 6'/2 > a'/2 and a',6' G (0,2). Note also that by Theorem 2.1 in 
[11], for all a;o = ao G [a, 6], almost surely 
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Therefore, the assumptions of assertion 1 of Corollary 3.2 hold with b{N) = {N + 

for any p' such that 1/p' G (0,1/6— l/min(2p,2)). And then, substituting in its proof 

Rn by i?^, almost surely 

sup sup N^/p'\R*M\<+oo, 

aG[o,6] 


which concludes the proof. 

C.2. Modulus of continuity and rate of convergence 

This section is devoted to the proofs of the results stated in Section 4.2. First, let us 
establish Theorem 4.2. 

Proof of Theorem 4.2. Let us fix xq = (oq, uq) G K^+i = [a, 6] x C (0,2) x 


Proof of assertion 1. Let us assume that p> b/2 and consider s the conditional pa¬ 
rameter defined by (6). Then, for any x = {a,u) G ATd-i-i and y = {a',v) G Kd+i^ 

s[x,y) <si[x,y) + S 2 {x,y), (38) 


where 


and 


'' +00 


1/2 


si(x,2/)= ^T-2/“|K(ai)-K(2/)|^ 


\n—l 


^ + CXD 


1/2 


S 2 (a:,y)= ^(T-V“-r-i/“') |K(y)|' 


First, let us focus on si. Note that for any x,y € A'd+i, 

si{x,y) < Cip{x,yf \og{l + p{x,y)~^)'', 

with Cl = (I]n“i Tn‘^^^\Ynf + ^ , where we have set 


Yn = sup 


mx)-Vn{y)\ 


p{x, yY log(l -f p{x, y)-YP ' 

x^y 


Since Kd is a convex compact set, applying a chaining argument and using the continuity 
of p, one checks that equation (13) implies that Yn G L^p. Then, since 2p> b>a and since 
the random variables Yn, n> 1, are i.i.d., Theorem 1.4.5 of [35] ensures that Ci < -foo 
almost surely. 
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Let us now focus on S 2 - Observe that \Vn{y)\ < with 

Xn = |14(a;o)| +CiYn 

for Cl = p{xo, zY\ log (1 + p{xci, . Let us remark that Ci < +cx), by con¬ 
tinuity of p on the compact set {xq} x Kd+i- Moreover, since Vnix^) € (X„)„>i is 

still a sequence of i.i.d. variables in LfP and following the same lines as for equation (36), 
we obtain that, almost surely, for any x,y & iLd-i-i, 

S2{x,y) < C2\a - a'\ 

with C 2 a finite positive random variable. Let us also note that by equation (8), there 
exist finite positive constants C2 and C3 such that for any x = (a, u) S Kd+i and any 
y={a',v) € Kd+i, 

|a - a'l < C2p{x,yf^^ < C3p{x,y) 

since H <1. Hence, since /3 G (0,1], almost surely, for any x,y G Xd+i, 
s(a:, y) < C'p(x, yf log(l -f p{x, 

with C a finite positive random variable. Then assertion 1 follows from Corollary 3.2. 

Proof of assertion 2. Let us choose p' > 0 such that 1/p' G (0,1/6— l/min(2,2p)). 
Then, replacing in the previous lines s by the parameter rjy and Theorem 1.4.5 of [35] 
by Theorem 2.2 of [11], we obtain: there exists C a finite positive random variable such 
that almost surely, for any x,y G Kd+i-, and for any N G'H, 

rNix,y) < p(a;,y)'^log(l -f 

Note also that by Theorem 2.1 in [11], almost surely 

sup ji?vr(xo)| < -foo. 

iVGN 

Therefore, by Corollary 3.4, almost surely, 

sup nGp' 

sup ji?Ar(a;)| <-foo, 

NeN xGKd+i 


which concludes the proof. □ 

Let us now prove Proposition 4.3. 

Proof of Proposition 4.3. Since equation (7) is fulfilled, there exists r G (0,1) such 
that p{x,y) < ho for all x,y G Kd+i with jjx — y\\ < r. Then, the assumptions done imply 
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that 


Xi := sup 

x,y&Kd+i 

0 <||a:-y||<r 


\Vi{x)-Vi{y)\ 
p{x,y)l^\log p{x,y)\'^ 


Gjh) 
h^iZho] m 


< sup 


:=G, 


where F{h) := h^\ \ogh\^. We assume without loss of generality that Hq = 2~^° with an 
integer fco > 1 is such that F is increasing on (0, ho] and equation (14) holds for h G (0, Hq]. 
Then, using the monotonicity of Q and F, 


q2p< 


+ 00 




(G{h)\ 




Therefore, by equation (14) and definition of F, 


+ 00 

E{X^P)<E{G^P)<max{2^,lfP ^ |fclog2|”^"^^'^ <+ 00 , 

k—ko 


which concludes the proof. □ 

C.3. Proof of Proposition 4.4 

Let Kd+i = [a,b] x 11^=1 C (0,2) x 

Proof of assertion 1. Let us fix an integer p > 1 and consider x^^l = {aj,ul^l) G 
Kd+i for each integer 1 < j < p. Then, we set x = (x^^l,..., x^pI). Choosing S = 
G Tn(^) > 0} we define Hg: (0, + 00 ) x <S x C ^ by 

Let us note that almost surely 

N 

Y, Hs{Tn,in,gn) = {.Sm,N . • ■ ■ , Sm{x^Pl)), 

n—1 

where Sm,N is defined by (4) with Wn given by (12). Then this series converges almost 
surely to (S'™ (cc^^i),..., 5^ {xIpI)). Since gi is symmetric, applying Theorem 2.4 of [34] 
and using a simple change of variables (t = rm{^)) and 12 (W^ \S) =0, we obtain that 

VA = (Ai,..., Ap) G CP, Vz e C, E(e'^(^^?=i = exp ils,x{z)), 

where 

IsA^) = I (^emnFMtXGg)))) 

J ( 0 ,+oo)xR‘* xC 
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- 1 - i5J(z(A, Jsit, (C,5)))l|D?(z(A.j5,(t.(c.g))»|<i)) dti/(dOFg(d9) 

with Pg the distribution of gi and 

Jsit, [f, g)) = f )5, • ■ •, ,f)g). 

Therefore, Is,\ does not depend on the function m, and then neither does the distribution 
of the vector Since this holds for any p and x, assertion 1 is 

established. 

Proof of assertion 2. Let us now consider the space B = C{Kd+i, C) of complex-valued 
continuous functions defined on the compact set Kd+i- This space is endowed with the 
topology of the uniform convergence, so that it is a Banach space. 

Let us assume that Sm belongs almost surely to 'Hp{Kd+i, (i^g) C B. For any x = 
{x^^\...,x^p1) G in view of its characteristic function, the vector ..., 

is infinitely divisible and its Levy measure is given by 


Fs{A)= / lA\{o}iHs{r,^,g))m{f)drv{d^)Fg{dg) 

J (0,+oo)x5xC 

for any Borel set A G B{Cp). We first assume that {a,u) i-)- fa{u,f) belongs to B for all 
^ G so that the function 

H : (0, -foo) X 5 X C —^ B, 

is well-defined. Since Hg is defined by (39), one checks that is a i3-valued 

infinitely divisible random variable with Levy measure defined by 


F{A)= lA\^o}{H{r,^,g))m{Odr,y{df)Fg{dg), A€B{B). 

J (0,+cxd) x»SxC 

Then, by Theorem 2.4 of [34], '^n^iF[{Tn,{fmgn)) converges almost surely in B as 
N —i> - 1 - 00 . Then, by definition of iL, the sequence (5'm,Ar)^gi^ converges in B almost 
surely. Therefore, its limit Sm is almost surely continuous on Kd+i- 

Let us now consider V C Kd+i a countable dense set in Kd+i- Then, since almost 
surely Sm G 'Hp{Kd+i, fi,g) and since Sm Sm, we get that almost surely 

\Smix) - Sm{y)\ 
p(x,2/)^[log(l -f p{x,y)-^)Y' 

x^y 


< -foo. 
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Then, by continuity of p, by almost sure continuity of Sm and by density of V on the 
compact set K^+i, 


|5'^(x) - Sm{y)\ 

x,yeX+i p(a;,y)/5[log(l + p(a;,y)-i)]'' °° 

x^y 

almost surely, that is, Sm belongs almost surely to Hpi^Kd+i, 13,7]). This establishes as¬ 
sertion 2 when {a,u) i—>■ is continuous for all ^ 

Assume now that (a, u) /a(u, 0 is continuous for €R‘^\Af with i^(A/") = 0 and set 

9a{u,0 ■= fa{u,O^R<‘\^r{0■ 

Then, almost surely, for all x = (a, u) S (0,2) x and all > 1, 

N 

Sm,Nix) = 

n—1 

and the conclusion follows from the previous lines since (a,u) ga{u,^) is continuous 

on Kd+i for all ^ The proof of Proposition 4.4 is then complete. 


Appendix D: Applications 

D.l. Proof of Proposition 5.1 


Let us first note that using Remark 5.1, we can and may assume without loss of generality 
that oi = 1, up to replace E by E/ai and te by 

Let us choose C > 0 arbitrarily small and consider the Borel function m defined on 

by 

w(C) = IICir“l||?||<A + r£;t(0“'’^^^|logTBt(0r^”‘’l||C||>A- 
Observe that m is positive on \ {0}. Then, 

0<c= / m(^)d^ = Cl-I-C 2 

with 


ci= [ ||Cir“d^ and C 2 =f rEt^O ''^^^|logri;t(^)| 


-i-C 


de 


Let us first observe that ci < oo since ao > 0. To prove that C 2 is also a finite constant, 
we need some tools given in [9, 31]. As in Chapter 6 of [31], let us consider the norm 
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II • list defined by 

lklUt=^ (40) 

According to the change of variables in polar coordinates (see [9]) there exists a finite 
positive Radon measure (Je* on Se^ = {^ £ : H^Ust = 1} such that for all measurable 

function if non-negative or in L^(]R‘^,d5), 

p p + OO p 

/ V^(C)d^= / / f{r^ 6)aEt{d9)r’^^^'^~^ dr. 

dRd Jo Js^t 

Applying this change of variables, it follows that C 2 < oo since C > 0. Hence, m = rh/c 
is well-defined and =m(^)df is a probability measure equivalent to the Lebesgue 

measure. Then we may consider 5'm(ao,u) defined by (15) for m £ so that Xa^ 
rfao5'm(ao, •) with dao given by (16). 

To study the sample path regularity of •) on Kd = 0^=1 apply Propo¬ 

sition 4.3 on Kd+i = {ao} x Kd C (0,2) x R'^ for 

Vi{ao,u) = 

with fata defined by (18). We recall that here is a random vector of R"^ with density 
771. Therefore let us now check that assumptions of Proposition 4.3 are fulfilled. 

For h> 0 and ^ £ R'^ we consider 

g{h,f)= min(cE* II , 1) l^ao (0 L 

where ce* > 0 is chosen such that |e'^“’^^ “ 1| < ce*\\te{u)^*^We*■ We consider the quasi¬ 
metric defined on R‘^+^ by 

p((a, u), {a', v)) = \a — a'\ pEiu,v), V(a,u), (a', u) £ R x R'^, 

which clearly satisfies equation (7). By definition of Vi, g and || • the random field 
Q = (g(h,fi))d^io -hoo) satisfies (i) and (ii) of Proposition 4.3. It remains to consider 
assumption (iii). Let 

I{h)=E{g{h)^)= [ (41) 

JRd. 

Since 'ipao satisfies (19), 

I{h)= [ mm{cEt\\h^"f.\\EtA) df < Ii{h)-i-hih) 

jR'i 

with 

Gm<A 
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where A is given by the condition (19), and 

Jr'^ 

From Lemma 3.2 of [7] there exists a finite constant Ci > 0 such that for all h G (0,e“^] 

Moreover, using again the change of variables in polar coordinates, there exists a finite 
constant C 2 > 0 such that for all h G (0,e“^], 

hih) < C2h2^|log(/i)|(^+^)("/“°-^^ 

Since /3 < oi, one find a finite constant Ca > 0 such that 

I{h) < (42) 

Hence, assumption (iii) of Proposition 4.3 is also fulfilled and applying this proposition, it 
follows that (13) is satisfied with j3 and 77 = l/ao + £j for all e > 0. Then, by Theorem 4.2, 
almost surely Sm G ’Hp(iFd+i, /5, l/Q;o + l/2 + e). By definition of p and this means 

that 


Sm{ctOi •) S 77p£ {Kd, /3,l/ao + l/2 + e). 

In particular, Smio^o,') is continuous on Kd- Then, since dagSmicto,') ^ag, Xag is 
stochastically continuous and almost surely 


C := sup 


l^agiu) >^Q|,(l')| 


u,veV,u^v Te{u — 7;)^[l0g(l + Te{u — t;)-l)]^/“‘’ + ^/^+® 


< +00, 


where V C Kd is a countable dense in Kd = 11^=1 foi- '^s write H* this event 

and let us define a modification of X^g on Kd- 

First, if w ^ H*, we set X^^{u){u}) = 0 for all u G Kd- Let us now fix w G H*. Then, we 
set 


Xl^{u){u)=Xg,g{u){w), 'iuGV. 

Let us now consider u G Kd- Then, there exists G T> such that lim„ ^+00 = U. It 

follows that, 

< C'(a;)rs(Mfo^ - w*^'”^)^[log(l+r£;(ufo^ - M^™^)”^)] ^ ^ 

so that (X*jj(ufoi)(w))^ is a Cauchy sequence and hence converges. We set 

Xl^{u){u:)= hin X:„(ufo))(cc). 

^ n—>-+cx) ^ 
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Remark that this limit does not depend on the choice of and that X*^(-)(a;) is 

then well-defined on K^. Observe also that, by stochastic continuity of Xag, X*^ is a 
modihcation of Xag. Moreover, by continuity of te. 


C{co) 


\X*^^{u){u;)-Xl^{v){co)\ 
u,v&Ki,u^v Te{u — 'u)/^[log(l -b Te{u — 


< -boo 


for all w G n and X* is continuous on K^. This concludes the proof. 


D.2. Multistable random fields 

This section is devoted to the proofs of the results stated in Section 5.2. Let us first 
establish Proposition 5.2. 

Proof of Proposition 5.2. Since p satisfies equation (7), so does p. Then, assumptions 
of Theorem 4.2 are fulfilled, which implies that {Sm,N) converges uniformly to Sm 
on Kd+i = [a, 6] x Kd- Therefore, iSm,N )converges uniformly to Sm on Kd since 
Sm,N{u) = Sm,N{o:{u),u) and Smiu) = Smi<^{u),u) and a is continuous. 

Moreover, by Theorem 4.2 there exists a finite positive random variable C such that 
for any u,v G Kd, 

|5m(u)<C'p(a;(M),a:(u))^[log(l-b/9(x(M),a:(u))"^)] ^ , 

where x(w) = (a(w),w). Moreover, by definition of p and since a € 7{p(Kd,l,0), there 
exists a finite positive constant ci such that 

\/u,vGKd, p{x{u),x{v)) <cip{u,v). 

Let us now recall that since p is continuous on the compact set Kd x Kd, M = 
sup„ p{u,v) < -boo. Then, up to change C, for all u,v G Kd, 

\Smiu) - Sm{v)\ <Cp{u,v)^ [log {1 piu,v)~^)] ^ 

since h i—log(l-b is increasing around 0 and bounded on [0,M]. 

Assertion 1 is then proved. Moreover, assertion 2 is a direct consequence of assertion 2 
of Theorem 4.2. The proof is then complete. □ 

Let us conclude this paper by the proof of Corollary 5.3. 

Proof of Corollary 5.3. Let Kd = 11^=1 [%: C and uq G Kd- Let us set 

a = mina, 6 = maxa and iLd+i = [a, 6] x C (0,2) x R'^. 

Kd Kd 
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Let us first note that Assumption 5 is fulfilled with Ki = [a, b] and then Sm is well-defined. 
Let us now consider pE and te as defined in Example 2.1. Then we set 

with C > 0 a parameter chosen arbitrarily small. Therefore, let us consider 

Vn{a,u) = fa{u,^n)'m{^n) 

where (^n)n>i is a sequence of i.i.d. random variables with common distribution /2(d^) = 
m(^)d^. The sequence {^n)n>i is assumed to be independent from (Tn,gn)n>i- Then, 
Assumption 4 is fulfilled. Moreover, 

\Vn{a,u) - Vn{a',v)\ < \Vn{a,u) - Vnia,v)\ + \Vnia,v) - Vnia',v)\. 


Let us set 


and 


Ci = 


\Vi{a,u) - Vi{a,v)\ 
Pe{u, v)\logpE{u, v)\'^ 

0 <||ii—?;||<r 


C 2 = sup sup 

uGKd 

ol^ol' 


\Vi{a,u) 


Vi{a',u)\ 


where r > 0 and the choice of 77 G R is given below. Then, for any x = (a, u) G Kd+i and 
any y = (o', v) G K^+i such that ||a; — y\\ < r. 


\V\{x) - Vi{y) \ < {Cl + C 2 ){pE{u,v)\\ogpE(u,v)\'' |q; - a'\) 

< Ci{Ci + C 2 )pix,y)\logp{x,y)\'^, 


where ci G (0,-|-oo) is a finite constant and p{x,y) = pe{u,v) + \a — a'\. Then, to apply 
assertion 1 of Proposition 5.2 with p = pe and /3 = 1, it suffices to establish that Ci,C 2 
and Vi{a{uo),uo) G (since b < 2). 

Let us first deal with Vi{a{uo),uo)- Using polar coordinates associated with E* (see 

[31]), 


/‘ + 00 2 

E(|Vi(a(uo),wo)| )<C 2 / min 11,1) di 

^0 


with C 2 G (0, -l-oo). Hence, Lemma 2.1 of [9] proves that Vi{a{uo),uo) G for any choice 

ofC 
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Let us now consider the random variable Ci. By homogeneity and continuity of tp, 
there exists a finite positive constant C 3 such that for any u,v € Kd, 

sup \Vi{a,u) - Vi{a,v)\ < - l|Zi 


with 

Zi = TEt (|i) ^ max (|log TEt iii )I, |log te* [fi) \ 

Combining the proofs of Propositions 4.3 and 5.1, we obtain that for any e > 0, choosing 
r small enough, 


E 




sup 


l|Zi 


1 u,v(iKd p£(M,v)|logp£;(M,w)|^^“^®. 
0<||it—i;||<r 


< + 00 . 


This implies that for any e > 0, Ci G for 77 = 1/a + e: and f well-chosen. 

Let us now study C 2 - Since Kd is a compact set, using polar coordinates and the Mean 
Value theorem, we have 

sup |Vi(a,7;) — bi(Q;',u)| < C 4 |a — a'\Z 2 


£/* 

with Z 2 = min(llr^jt ||, l) 2 ’i| logr£;t(^„) -|- C 5 I and C 4 and C 5 two finite positive con¬ 

stants. Using polar coordinates, one checks that Z 2 GLf, which implies that C 2 GLf . 

Therefore, for any e > 0, assumptions of assertion 1 of Proposition 5.2 are fulfilled for a 
well-chosen f. This implies that almost surely, for any e > 0, Sm G 'Hpj^{Kd, 1,1/a -I- 1/e) 
with a = minif_j a. Hence, for any e > 0, Sm G Kpj^^B(uo,r)(^o, 1, l/a(wo) + 1/2 -|- 1/e) for 
r small enough. This concludes the proof. □ 
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